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Using Quantum Monte Carlo we compute thermodynamics and spectra for the orbitally degenerate 
Hubbard model in infinite spatial dimensions. With increasing orbital degeneracy we find in the 
one-particle spectra: broader Hubbard bands (consistent with increased kinetic energy), a narrow- 
ing Mott gap, and increasing quasi-particle spectral weight. In opposition, Hund's rule exchange 
coupling decreases the critical on-site Coulomb energy for the Mott transition. The metallic regime 
resistivity for two-fold degeneracy is quadratic-in-temperature at low temperatures. 
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The Hubbard Hamiltonian, the simplest model for 
strongly interacting many-electron systems, has been ex- 
tremely popular in that, despite its simplicity, the model 
is considered to capture essential physics in electronic 
systems, ranging from a metal-insulator transition (MIT) 
and associated antiferromagnetism to possible d-wave su- 
perconductivity. Although most of the real systems dis- 
playing these phenomena have orbital degrees of free- 
dom, most theoretical works have concentrated on the 
orbitally non-degenerate model for simplicity. Recently, 
with the advent of the colossal magnetoresistance mate- 
rials 0], the proposed triplet pairing superconductivity 
in Sr2RuC>4 0], and alkali-doped fullerides (^^], atten- 
tion has been turned to the multi-orbital Hubbard model 
(MOHM) ||,|. 

With on-site orbital degeneracy iVdeg > 1, one has to 
consider two additional aspects of the problem: first, the 
effect of orbital degrees of freedom on the mobility of elec- 
trons and second, the role of on-site exchange interactions 
between degenerate orbitals. In this paper, we concen- 
trate on the MIT in MOHM motivated particularly by 
alkali-doped fullerides, A^Cqq (A=K, Rb, Cs etc.). We 
have studied the problem using Quantum Monte Carlo 
in a d — oo model on both Bethe and hypercubic lat- 
tices, with a focus on dynamical properties (obtained by 
analytic continuation of imaginary time data using the 
Maximum Entropy method [MEM]). We present the first 
systematic calculation of one-electron spectral functions 
at and away from particle-hole symmetry and the first 
calculation of the optical conductivity and d.c. resistiv- 
ity for the two-orbital model. 

We find that the upper and lower Hubbard bands dis- 
play a considerable broadening for N^eg = 2,3 relative 
to iVdog = 1, which supports the idea of Gunnarsson, 
Koch, and Martin jj| and Lu (6) that orbital degeneracy 
serves to increase the effective hopping matrix element 
and thereby mitigate the efficacy of correlations in driv- 
ing a MIT . We also find that the quasiparticle spectral 
weight increases with -/Vd eg . We find that the inclusion 
of Hund's rule exchange serves to decrease the effective 
hopping, thereby decreasing again the critical value of 
on-site Coulomb interaction U = U c needed for the MIT, 



according to the approximate expression 



U c (N doe , J) « y/NdegUcO.,0) ~ N deg J, 



(1) 



where J is the Hund's rule exchange energy. Finally, we 
present a calculation of the optical conductivity cr(u),T) 
and resistivity p(T) for the iVd eg = 2 case which gives a 
T 2 behaviour in p(T) for T -> 0. 

We were motivated to carry out this study by con- 
sideration of ^Ceo- Even though the local density ap- 
proximation predicted well the band-width and the band- 
gap of this system, it has been suggested that the strong 
Coulomb repulsion on the Ceo site could drive the sys- 
tem to a strongly correlated metal or even to a Mott- 
insulator 0. Since the band- width of A^G^q (0.4 eV) is 
more than two times smaller than the on-site Coulomb 
repulsion (1-1.5 eV), one can naively predict that the sys- 
tem should be a Mott-insulator, if interpreted in terms 
of an orbitally nondegenerate Hubbard model. However, 
the alkali-doped fullerides are metals despite the large 
resistivity at zero temperature jl0| . 

It has been suggested that the metallic A^G^q can 
be understood by introducing orbital degeneracy in the 
Hubbard model The MIT in the MOHM has been 

studied by the Gutzwiller approximation ||, where Lu 
showed that the critical value for the on-site Coulomb 
repulsion strength U c {N des ) for MIT at half-filling in- 
creases as (N dcg + l)U c with the orbital degeneracy iVdeg- 
This explains why alkali-doped fullerides (iVdeg = 3) can 
be metallic despite the large intra-molecular Coulomb 
repulsion and small conduction band- width jjj. Later, 
Gunnarsson et al. || argued that U c (N dcg ) behaves as 
\/Nd C g, in a much weaker fashion than the linear depen- 
dency of Ref. ||. The electron hopping is argued to be 
effectively enhanced with the orbital degeneracy due to 
the anti-ferromagnetic correlation between nearest sites, 
which is also supported by their numerical calculation. 

To carry out our studies, we have used the dynamical 
mean field theory (DMFT) jnj, where a lattice Hamilto- 
nian is rigorously mapped to an effective impurity model. 
This approximation, which is exact for infinite spatial di- 
mension or coordination number, is reasonable for mod- 



1 



elling the Ceo-derived t\ u band due to the large coordi- 
nation number, 12, between the Cgo molecules located at 
face-centered-cubic sites. Furthermore, with the on-site 
Coulomb interaction, the impurity Hamiltonian from the 
DMFT is expected to retain the essential physics of the 
MIT. The multi-orbital Hubbard Hamiltonian reads 

H = (tij4 m(T c jma + h.cj + e d n lma + H c , (2) 

ij,tna irtMT 

where t%j is the hopping matrix element between different 
sites i and j, m the orbital index, a the spin index, ed the 
d-level energy as the center of the non-interacting band, 
He the interaction Hamiltonian for Coulomb energy. For 
simplicity, we have neglected inter-band hopping since we 
concentrate on the local properties. 

The Coulomb interaction term is written (lj] as 

Hp = (U + J)/; n im ^n im i + (U - J) 2J n ima n im 'a 

m im<m'(T 
~\~U ^ ] n im'\ n im' i + J ^ ] C irn'^ C iml Cim ' l Cim T ■ (^) 

Here, we discard the last term proportional to J for com- 
putational convenience although, for full rotational sym- 
metry, all the terms in the above Hamiltonian are needed. 
The truncated Hamiltonian still captures the Hund's rule 
physics, and we believe that the form of Eq. (g) remains 
valid, possibly with a small correction to the term pro- 
portional to J. 

The effective impurity model in the DMFT is stud- 
ied by the Hirsch-Fye algorithm |l3| where the partition 
function is reformulated using a discrete path integral 
with an imaginary-time interval which we set to At = 
1/3. For the free density of states (DOS) in the DMFT, 
a semi-circular DOS ( corresponding to the infinite coor- 
dination number Bethe lattice), D(e) = §Vl e 2 with 
band-width W=2, is used for calculations of thermody- 
namic quantities such as occupation numbers. A Gaus- 
sian DOS (corresponding to the infinite dimensional hy- 
percubic lattice), D(e) — exp(— e 2 ) with W=y2, is 
used for one-particle spectral functions. We normalize 
tij = t by t = t*/2yd, as the dimensionality d goes to 
infinity. All energies here are measured in units of t* for 
the Gaussian DOS, and in units of the bandwidth for the 
semi-circular DOS. 

At J = 0, we identified the MIT for N dcg = 1,2,3 near 
half-filling by computing the occupation numbers at 
away from the particle- hole symmetric value, i.e., e d = 
e ph + 0.2 where e ph = {N deg - \)U + (1 - -§*) J. FIG. g 
shows the occupation numbers versus the Coulomb repul- 
sion with the semi-circular DOS. The occupation num- 
bers and the chemical potentials are plotted with respect 
to the particle-hole symmetric values. The temperatures 
are at T = 1/12,1/16. When the system is insulating 
with the chemical potential inside an insulating gap, the 
occupation number does not change with a slight shift 
of chemical potential from e p h- As will be illustrated 



later, the Mott insulating-gap manifests as a plateau in 
rid vs e<2 . The critical Coulomb repulsion, U c , is read 
off from the point where the occupation number rid at 
e<2 = e p h + 0.2 deviates from AT deg in FIG. |l|. 

Uc(Ndeg) increases with the orbital degeneracy, in 
agreement with Ref. ||, although the data are not suf- 
ficiently precise to conclude that the insulating gap is 
proportional to \/N deg from small N deg and U/W. For 
A^dcg = 3, the occupation numbers have not fully con- 
verged yet for T = 1/16 and we obtain U c /W > 2.3. 
Compared with A 3 C 60 , this value of U c puts the experi- 
mental ratio near the edge of MIT transition. Therefore, 
it is theoretically possible, if not decisively proven, that 
A3 Ceo can be described as a strongly correlated metal 
despite being at half filling. 

The QMC Green's function along imaginary time is 
analytically continued for the one-particle spectral func- 
tion by MEM One-particle spectral functions at 
and away from half- filling are shown in FIG. when 
T = 1/12 and J = 0. In plot (a), at half-filling, the 
upper (UHB) and lower Hubbard bands (LHB) are split 
by the Coulomb interaction strength U . Note that the 
width of these peaks increases with the orbital degen- 
eracy, in the same fashion that the electron hopping is 
enhanced for conduction of the charge excitation. For the 
non-degenerate model (iVdcg = 1), the spectral weight at 
zero frequency is very small, hence a sign of an insulating 
pseudo-gap (the system cannot become a true insulator 
due to the non- vanishing tail of the Gaussian free DOS). 
With increasing Ad cg , the MOHM evolves to solutions 
with more metallic character. 

The peak positions of the UHB and LHB are progres- 
sively shifted further away from ±Z7/2 as iVdeg increases. 
These shifts, roughly proportional to N dcE t 2 /U for U suf- 
ficiently large, can be understood in terms of the pertur- 
bational contributions to energy levels with occupation 
numbers iVdeg - 1 5 Ad cg and -/Vdeg + 1, respectively. When 
t w U , one should also take into account this shift in 
addition to the band broadening at the MIT, which can 
complicate the argument of Gunnarsson, Koch, and Mar- 
tin |5j that led to a ^Ndcg dependence to U c . 

The central peaks for 7Vd og = 2,3 are the quasi-particle 
(QP) excitation which emerges from the scattering of the 
spin degrees of freedom of conduction electrons. The 
width of the QP peak defines the new low temperature 
energy scale (renormalized Fermi energy). Since the spec- 
tral weight from the UHB and LHB is expected to be 
larger at ui = for larger iVd eg due to the broadening 
mentioned in a previous paragraph, the QP peak accord- 
ingly has larger weight at the chemical potential. This 
trend goes with enhanced metalicity for increasing N deg . 
We will discuss more the dependency of the QP weight 
upon orbital degeneracy in relation to the electronic re- 
sistivity. 

FIG. |^ (b) shows the one-particle spectral functions 
away from half-filling at U = 4, ed = — 1 for all iVd cg . In 
the paramagnetic phase, the total occupation number is 
evenly distributed among the orbitals and the spectral 
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weight of the LHB roughly goes as 1/N dcg . Since \ed\ <C 
U and the doubly occupied configuration exists in small 
quantum weight, total occupation number is not very 
sensitive to the orbital degeneracy. 

The behavior of the MIT with a finite Hund's rule cou- 
pling J is shown in FIG. ||. When the exchange interac- 
tion J in Eq. (^) is included, the Coulomb repulsion is 
effectively increased, favoring the insulating phase. Nu- 
merical calculation confirms that the insulating phase 
is expanded with finite positive J. Figure || shows the 
plot of n d vs. e d — e ph with J = 0.0, 0.2, 0.5, 1.0 at 
Adcg = 2 and = 16. The total occupation num- 
ber rid begins to deviate from the half-filling value A^ eg 
when the system become metallic, e^'s where the metal- 
lic phases begin are pushed to higher values with Ae d 
approximately linear in J. In other words, the insulating 
gap has increased with finite J. 

The phenomena can be easily understood by consider- 
ing the change in the on-site energy at half-filling. The ef- 
fective on-site Coulomb repulsion U e f / can be computed 
as E(N dcg + 1) + E(N dcg - 1) - 2E(N dcg ) from the on-site 
Hamiltonian Eq. (||), which leads to 

U ef f = U(N dcg , J = 0) + N dcg J, linear in N dcg . (4) 

Also, we have to remember that the hopping of elec- 
trons leaves behind a trail of broken Hund's rule mul- 
tiplets. Therefore, there would be additional effective re- 
pulsion due to hopping, which is ignored in Eq. (Q). For 
^Ceo, J is one order of magnitude smaller than U 0, 
the Hund's rule coupling plays less important role than 
^/N dcg U c (l, 0) in Eq. (||). For systems with larger J, 
one should take into account the last term in Eq. (|3|) to 
properly describe the magnetic correlation. 

Finally, we discuss the transport properties of the 
MOHM. From the analytically continued one-particle 
spectral function, we can extract the electron scatter- 
ing rate (imaginary part of the self-energy) and thereby 
calculate the optical conductivity, er(oj), within the 
DMFT Jl|l, as shown in FIG. f§ The inset (a) shows the 
evolution of spectral functions for a half-filled doubly de- 
generate model at inverse temperatures (3 = 2, 4, 8, 12, 16, 
where the curve with the large QP peak (labelled as a in 
the figure) corresponds to low T. The optical conduc- 
tivities shown in the main plot are from the same set 
of temperatures. The optical conductivity is computed 
by convoluting the one-particle spectral functions and 
is therefore characterized by transitions between various 
single-particle excitations, as indicated in the figure. The 
Drude peak (marked as a — a), resulting from the scat- 
tering within the QP peak a, diverges as the temperature 
goes to zero. Note that, since the leading order of <t(uj) 
is 0(l/d), we have plotted da (us). 

The resistivity p(T), plotted in inset (b), is calcu- 
lated from the ix> — * limit of the optical conductivity, 
piT)- 1 = lim^-,0 cr(w). The unit of p(T)/d is pflcm @. 
The data shown here are for half-filled Hubbard models 
with double (filled circle) and triple (open circle) orbital 



degeneracy at U = 4 using the Gaussian free DOS with 
effective band width W = \pl. The resistivity for doubly 
degenerate model shows a cross-over to a Fermi liquid 
behavior, p{T) oc T 2 , at low T from a resistivity of in- 
sulating character at high T. For the triply degenerate 
model, we could not determine the low temperature law 
of resistivity due to noise in QMC data. 

It should be noted that p(T) for iVdeg = 3 is about 
one order of magnitude smaller than p(T) for N dcg = 2, 
in contrast to the weak iVd C g-dependency in the low-L 7 
perturbation theory |17|. Although not discussed in this 
paper, the magnitude of p(T) seems to strongly depend 
upon the shape of band edge of free DOS. For exam- 
ple, the resistivity for semi-circular DOS is order(s) of 
magnitude larger than for Gaussian DOS at the same 
U/W ratio. A systematic investigation of resistivity on 
the orbital degeneracy and DOS and its comparison with 
experiments are left for further studies. 
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FIG. 2. Spectral lunctions for N dcg = 1, 2, 3 at U = 4 and 
J = 0. (a) The effective band-width of the charge excitation 
peaks (nearly at ±U /2 at half-filling) is broadened with in- 
creasing iV dc g due to the enhanced hopping for one-particle 
(hole) doped states. The quasi-particle (QP) peak becomes 
dominant for large AT deg , as more weight of charge excitation 
peaks is transferred to the QP peak, (b) At a fixed t d = — 1, 
all three cases have total occupations near 0.9 and the inte- 
grated weight up to the chemical potential is roughly 1/N deg . 
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FIG. 1. The occupancy rid versus the Coulomb repulsion 
U at fixed = e p h — (Ndcg — §)t/+0.2. The critical value of U 
(marked with arrows at T = 1/16) at the cross-over between 
insulating and metallic solutions increased with iVdeg- The 
temperature, T = 1/16, is not low enough for iVdcg = 3 case to 
equate the U c = 2.3VK as the zero-temperature quantum crit- 
ical point. Due to the weak functional dependency of y* 1 N dcg 

and small A?dcg, it is not clear that U c (Ndeg) = \J -/Vdcg £4(1). 
(See inset.) Lines are guide to the eye. 
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FIG. 3. n d vs. e d — e vh at J = 0.0, 0.2, 0.5. For positive J 
the onset of a metallic solution is shifted to larger values of td 
from the particle-hole symmetric parameter, e p h- The shift is 
roughly linear in J, leading to a suggestion that the on-site 
Coulomb interaction is responsible with AU(J) = NdcgJ- 
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FIG. 4. Optical conductivity and resistivity for Ndeg = 2 
at U = 4, e<j = — 1. The optical conductivity is made up of 
three excitations arising from transitions between three struc- 
tures in the spectral functions plotted in the inset (a). The 
spectral functions at /3 = 2,4,8, 12, 16 have QP peak, LHB, 
UHB marked as a, b, c, respectively. The Drude peak, arising 
from the scattering within the QP peak, becomes stronger 
as T lowers. The d.c. resistivity p(T) for the doubly degen- 
erate-band Hubbard model (filled circle) illustrates that the 
Fermi liquid behavior (dashed line) holds at half-filling below 
the critical U value. The triply degenerate case (empty circle) 
has much lower values for the resistivity. The unit of p(T) /d 
in the figure is (iQcm |16j. 



5 



